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Search for Gluinos and Scalar Quarks in pp Collisions at ^/s = 1.8 TeV using the 

Missing Energy plus Multijets Signature 
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(Dated: February 7, 2008) 
We have performed a search for gluinos (g) and squarks (q) in a data sample of 84 pb~^ 
of pp collisions at -^s = 1.8 TeV, recorded by the Collider Detector at Fermilab, by 
investigating the final state of large missing transverse energy and 3 or more jets, a 
characteristic signature in _R-parity-conserving supersymmetric models. The analysis has 
been performed 'blind', in that the inspection of the signal region is made only after the 
predictions from Standard Model backgrounds have been calculated. Comparing the data 
with predictions of constrained supersymmetric models, we exclude gluino masses below 
195 GeV/c^ (95% C.L.), independent of the squark mass. For the case niq ~ mg, gluino 
masses below 300 GeV/c^ are excluded. 
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The Standard Model(SM) @] accurately describes 
physical phenomena down to scales of ~ 10""'^^ cm. 
There are many extensions of the Standard Model to 
smaller length scales, including extra gauge interac- 
tions, new matter, new levels of compositeness, and 
supersymmetry (SUSY). Of these, supersymmetry |^ 
treats the bosonic and fermionic degrees of freedom 
equally and provides a robust extension to the Stan- 
dard Model . For simplicity the minimal construction 
(MSSM) is often used to link SUSY with the Stan- 
dard Model The most general MSSM would in- 
duce proton decay with a weak-interaction lifetime; to 
avoid this, baryon and lepton conservation are enforced 
in the MSSM by postulating a new conserved quan- 
tity, i?-parity, R — (_i)3(s-L)+2s^ where for each par- 
ticle s is the spin, and B and L are the respective 



baryon and lepton assignments. i?-parity conservation 
leads to characteristic SUSY signatures with missing 
transverse energy in the final state due to the stable 
lightest supersymmetric particle (LSP). We assume in 
the search described below for the bosonic partners of 
quarks (squarks) and the fermionic partners of gluons 
(gluinos) that the LSP is weakly interacting, as is the 
case for most of the MSSM parameter space. 

We consider gluino and squark production within the 
minimal supergravity model (mSUGRA)]^. In this 
model the entire SUSY mass spectrum is essentially 
determined by only five unknown parameters: the com- 
mon scalar mass at the GUT scale, Mq; the common 
gaugino mass at the GUT scale, M1/2] the common 
trilinear coupling at the GUT scale, Aq; the sign of 
the Higgsino mixing parameter, sign(fj,); and the ratio 
of the Higgs vacuum expectation values, tan/?. Min- 
imal SUGRA does not make predictions for the part 
of the rrig-mg mass parameter space where squarks of 
the first two families are lighter than about 0.8 times 
the mass of the gluino. Hence for ruq < m,g we use the 
constrained MSSM ||^ with the set of input parameters 
being the mass of the gluino, m^; the CP-odd neutral 
scalar Higgs mass, niA] the squark masses, Tig^; the 
slepton masses, ; the squark and slepton mixing 
parameters, ; and and tan/3. 

We investigate whether the production and decay of 
gluinos and scalar quarks is observable in the rate of 
>3-jet events with large missing transverse energy at 
the Collider Detector at Fermilab (CDF). The large 
missing energy would originate from the two LSPs 
in the final states of the squark and gluino decays. 
The three or more hadronic jets would result from 
the hadronic decays of the q and/or g. We use the 
ISAJET Monte Carlo (MC) program Q with tan /^ = 3 
to generate datasets of squark and gluino events, and 
the PROSPINO program ^ to calculate the production 
cross sections. To be conservative, only the first two 
generations of squarks [u, d, c, s) are assumed to be pro- 
duced |6) in the general MSSM framework; we addition- 
ally consider production of the bottom squark (6) in the 
mSUGRA case. The search is based on 84 ± 4 pb~^ of 
integrated luminosity recorded with the CDF detector 
during the 1994-95 Tevatron run. 

The CDF detector is described in detail else- 
where 0] . The momenta of charged particles are mea- 
sured in the central tracking chamber (CTC), which is 
positioned inside a 1.4 T superconducting solenoidal 
magnet. Outside the magnet, electromagnetic and 
hadronic calorimeters arranged in a projective tower 
geometry cover the pseudorapidity region \ri\ < 4.2 ^ 
and are used to identify jets. Jets are defined as lo- 
calized energy depositions in the calorimeters and are 
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FIG. 1: The ^ spectrum after the online trigger and 
the two stages of the data preselection. The numbers of 
events surviving the first and second selections are 892,395 
and 286,728, respectively. The variables Equt, Nqut are 
energy and number of towers out of time [Q. 

reconstructed using an iterative clustering algorithm 
with a fixed cone of radius AR = \/ Ar\^ + Z\^2 = q.T 
in 77— space . Jets are ordered in transverse energy, 
Et — E sin 9, where E is the scalar sum of energy de- 
posited in the calorimeter towers within the cone, and 
6 is the angle formed by the beam-line, the event ver- 
tex and the cone center. 

The missing transverse energy is defined as 
the negative vector sum of the transverse energy 
in the electromagnetic and hadronic calorimeters, 
= — sin0i)ni, where Ei is the energy of the 

«-th tower, fii is a transverse unit vector pointing to the 
center of each tower, and 9i is the polar angle of the 
tower; the sum extends to |ry| < 3.6. The data sam- 
ple was selected with an on-line trigger which requires 
li- = llrl > 30 GeV. 

We use a two-stage preselection to reject accelerator- 
and detector-related backgrounds, beam halo, and cos- 
mic ray events. The first stage is based on timing 
and energy information in the calorimeter towers to re- 
ject events out-of-time with a pp collision. The second 
stage uses the event electromagnetic fraction (Fem) and 
event charged fraction {Fch) to distinguish between real 
and fake jet events The preselection requirements 
and the corresponding missing transverse energy spec- 
tra are presented in Figure 1. At least three jets with 
Et > 15 GeV, at least one of them within jr/j < 1.1, 



are then required in events that pass the preselection. 
A total of 107,509 events, predominantly from QCD 
multijet production, survive the three-jet requirement. 

The observed missing energy in QCD jet production 
is largely a result of jet mismeasurements and detector 
resolution. In a QCD multijet event with large miss- 
ing energy, the highest Et jet is typically the most 
accurately measured. When the second or third jet is 
mismeasured because it lands partially in an uninstru- 
mented region (a 'gap'), the fjT is pulled close in (j) to 
the mismeasured jet. A jet is considered non-fiducial 
if it is within 0.5 rad in of the fr direction and also 
points in to a detector gap. The second and third 
highest Et jets in an event are required to be fiducial. 
We eliminate the residual QCD component by using 
the correlation in the S4>i = | (pleading jet — '/'^l ver- 
sus (502 = l^socond jet " 4'^ \ plauc. We accept events 
with Ri > 0.75 rad and R2 > 0.5 rad, where Ri = 
^S(bl + {n-6(bi)^and R2 = ^ b<t>\ + - 502)2. 

To avoid potential a 'posteriori biases when searching 
for new physics in the tails of the missing transverse 
energy distribution, once we define the signal candi- 
date data sample we make it inaccessible. This anal- 
ysis approach is often referred to as a 'blind analysis' 
and the signal candidate data sample as a 'blind box'. 
The 'blind box' data are inspected only after the entire 
search path has been defined by estimating the total 
Standard Model backgrounds and optimizing the sen- 
sitivity to the supersymmetric signal. We use three 
variables to define the signal candidate region : fir-, 
Ht = Et{2) + Et{z) + Pr, and isolated track multi- 
plicity, Nl^f. [|lj]. The 'blind box' contains events with 
It > 70 GeV, Ht > 150 GeV, and Nl^^=0. The large 
missing transverse energy requirement for the defini- 
tion of the box is motivated by the requirement that 
the trigger be fully efficient The Ht requirement 
provides good discrimination between signal and back- 
ground ||l^. The Nl^^ requirement increases the sen- 
sitivity of the search for all-hadronic final states by 
significantly reducing the backgrounds from W/Z+iets 
and top-antitop (tt) events while retaining the signal 
cascade decays in which a lepton is produced close to a 
jet (non- isolated lepton). The analysis path is shown in 
Table |. We reduce the background contribution from 
W{-^ ej^)+jets and tt production by requiring the two 
highest energy jets not be purely electromagnetic (jet 
electromagnetic fraction f^m <0.9). We further reduce 
the contribution from QCD backgrounds (mismeasured 
jets) by requiring the fr vector not be closer than 0.3 
rad in to any jet in the event. 

We estimate the W and Z boson backgrounds 
by using a leading order perturbative QCD calcu- 
lation for W(Z)+ jets as implemented in the VEC- 
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TABLE I: The data selection path for the fh- +> 3 jets 
search. After the fourth step, all events that could fall in the 
'blind box' are removed from the accounting. The events 
tabulated in the following steps are only in the control bins. 



Requirement 


Events 


Preselection 


286,728 


Njet >3 {AR = 0.7, Et> 15 GeV) 107,509 


Fiducial 2nd, 3rd jet 


57,011 


Ri >0.75 rad,i?2 >0.5 rad 


23,381 


|i^>70 GeV, Ht>150 GeV, 






'blind box' 


Et(i) >70GeV 




Et(2) >30GeV 




|77|(1 or 2 or 3) < 1.1 


6435 


fem(l) , /em(2) < 0.9 


6013 


L2 trigger requirement 


4679 


5(t)min{$T — jet) > 0.3 rad 


2737 



BOS Monte Carlo |T^, enhanced with a coherent par- 
ton shower evolution of both initial- and final-state 
partons, hadronization, and a soft underlying event 
model (VECBOS-I-HERWIG @). Events with large 
missing transverse energy and >3 jets in the final state 
are expected primarily from Z{—f vi>)-\- >3 jets and 
W{-^ Tv)+ >2 jets (the third jet originating from the 
hadronic r decay) processes. The MC predictions for 
events with > 3 jets are normalized to the observed 
Z{-^ ee)-\- jets data sample via the measured 



ratio, where Nj^t is the number of jets. The ratio 
p = Tl^^^T^~r'^-t'f^-*t^ ^ is used to normalize the W 

^ (T{pp^Z{^e+ e )-\-jets) 

MC predictions. Assuming lepton universality, the pre- 
dictions for the number of events with > 2- and > 3- 
jets from W and Z production and decay to all fla- 
vors are normalized to the data for Z(— > e+e~)-|- > 2 
jets. By normalizing the MC predictions to data we 
avoid large systematic effects due to the renormal- 
ization scale, the choice of parton density functions, 
initial- and final-state radiation, and the jet energy 
scale. The total uncertainty (~10%) is then dominated 
by the uncertainty on the luminosity measurement, the 
uncertainty on the measured ratio jv^^t+i ' ^'^"^ 
certainty on the predicted ratio p as a function of Njet- 
We estimate the backgrounds from single top, tt, 
and diboson events with Monte Carlo predictions nor- 
malized using the respective theoretical cross sec- 
tion calculations for these processes. We generate 
tt events with the pythia MC program nor- 
malizing to the fully resummed theoretical cross sec- 
tion a^^ = 5.06i|^:3^ Pb for rutop = 175 GeV/c^ 
p8[ |. We assign a total uncertainty of ±18% on the 
cross-section to take into account the uncertainty on 
the top quark mass. The top quark can also be 



TABLE IL Comparison of the Standard Model prediction 
and the data in the bins neighboring bin 8, the 'blind box'. 
After the contents of cons were compared in detail to stan- 
dard model predictions, we 'opened the box'. We find 74 
events in bin 8. 



Bin Definition 


EWK 


QCD 


All 


Data 


fh- >70,HT>150,Nl^°k> 


14 


6.3 


20±5 


10 


fJr >70,JJT<150,iVtV'fc= 


2.3 


6.3 


8.6±4.5 


12 


35<^ <70,i?T>150,AftV'fc= 


1.95 


135 


137±28 


134 


fJr >70,i?T<150,iVt"fc> 


1.73 


<0.1 


1.73±0.3 


2 


35<^ <70,HT>l50,Nlrk> 


14 


9.4 


23.4±6 


24 


35<|t <70,HT<l50,Nlrk^0 


5 


413 


418±69 


410 


35<|t <70,HT<150,Nli'°t:> 


3.3 


28 


31±10 


35 


>70,//T>150,iVtV'fc= 


35 


41 


76±13 


n 



produced singly via W^-gluon fusion and qq annihi- 
lation with cross sections of awg = 1-7 pb (± 17%), 
and ayy.^ti = 0.73 pb(± 9%) [|T8). We use the her- 
WIG (VF-gluon fusion) and pythia {qq annihi- 
lation) programs to generate the single top produc- 
tion processes. We generate boson pair production 
with the PYTHIA MC and use the calculated cross 
sections aww = 9-5 ± 0.7 pb, awz — 2.6 ± 0.3 pb and 
azz = 1.0 ±0.2 pb ll|. 

The data samples we use to study and normalize the 
QCD Monte Carlo predictions consist of events col- 
lected by on-line identification of at least one jet with 
transverse energy above trigger thresholds of 20 and 
50 GeV, and with integrated luminosity of 0.094 pb~^ 
and 2.35 pb~^, respectively. The corresponding QCD 
MC samples are generated using the HERWIG program 
and a CDF detector simulation. The shapes of the frp 
and jet multiplicity distributions are in good agreement 
with the data, as are the jet kinematic distributions. 
The QCD predictions are absolutely normalized to the 
data for Nj^t > 3. The total uncertainty on the QCD 
background estimate is ^15%, dominated by a 12% 
uncertainty due to the detector resolution. 

There are seven bins around the 'blind box' formed 
by inverting the requirements which define it {i.e. by 
changing the direction of the inequalities shown in the 
bin definitions of Table ||). We compare the Standard 
Model background predictions in the bins around the 
'blind box' with the data. The results are shown in 
Table Of the 35 events from electroweak processes 
predicted in the 'blind box', ^^37% are expected from 
Z vV+ > 3 jets, -20% from W tv+ > 2 jets, 
^^20% from the combined W e(/i)i/e(i'^)+ > 3 jets, 
and ^--^20% from tt production and decays. We also 
compare the kinematic properties between Standard 
Model predictions and the data around the box and 
find them to be in agreement [p^ . 

To probe the SUSY parameter space in a simple and 
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FIG. 2: Comparison in the 'blind box' between data 
(points) and Standard Model predictions (histogram) of 
fJr, Njet, leading jet Et and Ht distributions. There are 
74 events in each of these plots, to be compared with 76±13 
SM predicted events. Note that the f}r distribution is plot- 
ted with a variable bin size; the bin contents are normalized 
as labelled. 

comprehensive way we divide the rriq — nig plane into 
four general regions : (A) > trig (mSUGRA, five 
q); (B) rriq rrig (mSUGRA, five q); (C) < m-g 
(MSSM, four q); (D) < rUg (MSSM, four q). We 
analyze representative points of each region and opti- 
mize the frr and Ht requirements for increased sensi- 
tivity to the signal using MC data. The ratio "^^^^ 
is maximized in region A for > 90 GeV and Ht> 
160 GeV; in region B for > 110 and Ht> 230 GeV; 
in C for > 110 and Ht> 170 GeV; and in D for 
> 90 and Ht> 160 GeV, where Nsusy is the num- 
ber of signal events and Nsm is the number of Stan- 
dard Model background events. The signal efficiency 
ranges between 1% and 14% for the different points in 
the parameter space, and its total relative systematic 
uncertainty (mostly due to parton density functions, 
gluon radiation, renormalization scale and jet energy 
scale) ranges between 10% and 15%. 

In the 'blind box', where we expect 76±13 Standard 
Model events, we observe 74 events. In Figure 2 the 
predicted Standard Model kinematic distributions are 
compared with the distributions we observe in the data. 
For the A/D, B and C region requirements, we ob- 
serve 31, 5 and 14 events where we expect 33 ± 7, 3.7 
± 0.5 and 10.6 ± 0.9 events respectively. Based on the 



FIG. 3: The 95% C.L. limit curve in the — mj plane 
for tan/3 — 3; the hatched area is newly excluded by this 
analysis. Results from some previous searches are also 
shown(CDF II, D0 |2|], LEP I [|^; the area at lower 
masses in the plane has been previously excluded by the 
UAl and UA2 experiments | p5| , p^ . The region labelled 
as rriq < m~o is theoretically forbidden as the squarks are 
predicted to be lighter than the LSP. 

observations, the Standard Model estimates and their 
uncertainties, and the relative total systematic uncer- 
tainty on the signal efhciency, we derive the 95% C.L. 
Jl9[ upper limit on the number of signal events. The 
bound is shown on the niq — rrig plane in Figure 3. For 
the signal points generated with mSUGRA the limit is 
also interpreted in the Mq — A/1/2 plane jl^. Studies 
of the dependence on the value of tan /? can be found 




In conclusion, a search for gluinos and squarks 
in events with large missing energy plus multi- 
jets excludes at 95% C.L. gluino masses below 300 
GeV/c^ for the case rriq w rrig, and below 195 GeV/c^, 
independent of the squark mass, in constrained super- 
symmetric models. This is a significant extension of 
previous bounds. 
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